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ABSTRACT 

The  temperature  driven  gas  flow  in  a  two-dimensional  fi¬ 
nite  length  MicroChannel  and  a  cylindrical  tube  are  studied 
numerically  with  the  goal  of  performance  optimization  of  a 
nanomembrane-based  Knudsen  compressor.  The  numerical  so¬ 
lutions  are  obtained  using  direct  simulation  Monte  Carlo  method 
and  discrete  ordinate  method  for  BGK  model  kinetic  equation  in 
a  wide  range  of  Knudsen  numbers  from  0.05  to  50.  The  length- 
to-height  ratios  from  5  to  30  were  examined.  Three  different  wall 
temperature  distributions  were  considered,  namely,  linear,  step¬ 
wise,  and  a  non-monotonic  profile  typical  for  a  radiantly  heated 
Knudsen  compressor  membrane.  The  short  channel  end  effects 
are  characterized,  and  the  sensitivity  of  the  mass  flow  rate  to  a 
non-monotonic  temperature  distribution  is  shown. 

INTRODUCTION 

Lab  on  a  chip  technology  is  drawing  the  attention  of  scien¬ 
tists  from  many  disciplines.  Recent  advances  in  MEMS  man¬ 
ufacturing  have  made  it  possible  to  construct  microscale  ana¬ 
lytical  sensors  such  as  integrated  gas  chromatography  systems, 
miniature  spectrometers,  and  mass  spectrometers.  The  miniatur¬ 
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ized  detection  devices  will  require  microscale  roughing  pumps  to 
provide  sensor  elements  with  gas  samples  at  necessary  environ¬ 
mental  conditions.  A  pumping  mechanism  that  can  be  exploited 
at  microscale  is  thermal  transpiration,  a  rarefied  gas  effect  that 
drives  the  gas  flow  along  the  temperature  gradient  in  a  tube  or 
channel.  The  main  goal  of  this  paper  is  the  numerical  study 
of  thermal  transpiration  flows  in  short  microchannels  to  aid  in 
the  performance  optimization  of  a  transpiration  based  microscale 
roughing  pump  known  as  the  Knudsen  compressor  [1], 

In  early  1900s  M.  Knudsen  built  and  studied  the  first  transpi¬ 
ration  based  compressor,  consisting  of  a  series  of  differentially 
heated  and  cooled  capillaries.  Each  stage  of  the  compressor  had 
a  capillary  section  where  the  wall  temperature  increases,  caus¬ 
ing  a  thermomolecular  pressure  build-up  at  the  high-temperature 
end  of  the  capillary.  The  capillaries  are  followed  by  a  connec¬ 
tor  stage  of  a  significantly  larger  cross-sectional  area,  where  the 
pressure  is  almost  constant  and  temperature  is  decreased  to  its 
original  value  at  the  beginning  of  the  stage.  The  modern  ver¬ 
sion  of  Knudsen  compressor  was  suggested  by  Pham-Van-Diep 
el  al  [2]  and  demonstrated  by  Vargo  and  Muntz  [3].  The  most 
critical  element  of  the  Knudsen  compressor  developed  at  USC 
is  the  thermal  transpiration  membrane  made  of  porous  materials, 
such  as  aerogel,  with  pore  diameters  on  the  order  of  the  mean  free 
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path  of  the  gas  at  standard  atmospheric  conditions.  The  temper¬ 
ature  gradient  is  maintained  across  the  transpiration  membrane 
by  resistive  or  radiant  heating.  Recently,  a  single-chip  micro- 
machined  implementation  of  a  Knudsen  pump  was  reported  in 
Ref.  [4],  which  can  evacuate  a  cavity  to  0.46  atm  while  operating 
at  atmospheric  pressure  and  using  80mW  input  power. 

Knudsen  compressor  has  several  attractive  features,  such  as 
no  moving  parts  and  no  fluids  or  lubricants,  and  therefore,  can 
have  a  number  of  promising  applications  provided  that  its  per¬ 
formance  can  be  optimized  to  meet  the  existing  pumping  re¬ 
quirements.  The  two  most  important  performance  parameters  of 
the  compressor  are  the  energy  use  and  volume  per  unit  through¬ 
put.  To  accomplish  minimization  of  these  performance  parame¬ 
ters  the  Knudsen  compressor  must  also  be  designed  to  maximize 
the  gas  conductance  through  the  device  for  a  given  temperature 
change.  The  required  flow  rate  in  the  Knudsen  compressor  can  be 
obtained  by  adjusting  the  cross-sectional  area  of  the  stages  and 
the  temperature  distribution  along  the  transpiration  membrane  at 
different  stages.  The  main  element  of  the  transpiration  based 
pump  is  rarefied  gas  flow  through  a  capillary. 

Rarefied  gas  flows  of  different  gases  through  long  tubes  and 
channels  were  a  topic  of  an  extensive  experimental  and  theoret¬ 
ical  research  over  several  decades,  and  a  comprehensive  review 
of  the  numerical  data  and  analytical  results  for  rarefied  gas  flows 
in  capillaries  was  performed  by  Sharipov  and  Seleznev  [5]  (see 
also  references  therein).  There  is  however  a  lack  of  reference 
data  for  nonisothermal  flow  through  a  capillary  of  a  finite  length, 
especially  for  arbitrary  pressure  and  temperature  drops. 

The  present  numerical  modeling  of  thermal  transpiration 
flows  in  microchannels  addresses  several  questions  important  for 
design  optimization  of  Knudsen  compressors,  namely:  the  in¬ 
fluence  of  different  temperature  internal  distributions  across  the 
membrane  on  the  flow  field  inside  the  capillaries  and  the  total 
transpiration  mass  flow  rates,  the  impact  of  the  rarefaction  on  the 
capillary  throughput,  and  capillary  end  effects. 

Problem  Statement 

In  the  present  work,  two  reservoirs  filled  with  a  single¬ 
species  gas  and  joined  by  a  two-dimensional  capillary  or  a  cylin¬ 
drical  tube  of  a  finite  length  L  and  height  h  or  radius  a  are  consid¬ 
ered  in  the  computations.  The  two-dimensional  channel  config¬ 
uration  corresponds  to  the  case  of  a  capillary  of  infinite  width  w. 
A  schematic  of  the  geometric  configuration  and  notations  used 
throughout  this  paper  are  shown  in  Fig.  1.  The  subscripts  I  and 
II  refer  to  the  quantities  attributed  to  the  left  and  right  container, 
respectively  (far  enough  from  the  channel  inlet  and  outlet).  The 
channel  walls  have  some  temperature  distribution  varying  from 
7/  to  7//.  If  the  reservoirs  are  at  the  same  temperature  (Tj  =  77/) 
and  the  pressure  pj  is  larger  than  pu  there  will  be  a  mass  flow 
from  the  left  container  to  the  right  one.  Conversely,  if  the  pres¬ 
sures  are  equal  but  a  temperature  drop  is  maintained  between  the 


reservoirs  then  the  gas  flow  develops  from  the  cold  container  to 
the  hot  one. 

The  summary  of  cases  that  have  been  considered  is  given  in 
Table  1.  The  gas  is  molecular  nitrogen  in  all  simulations.  The 
first  set  of  computations  (Cases  I-a  to  I-e)  is  for  a  flow  with  both 
pressure  and  temperature  gradients  in  a  channel  with  a  length- 
to-height  ratio  of  5  at  different  Knudsen  numbers.  Hereafter,  the 
Knudsen  number  is  based  on  the  mean  free  path  in  the  right  reser¬ 
voir  and  the  channel  height  or  tube  diameter.  Note  also  that  the 
change  in  Knudsen  number  is  due  a  change  in  the  mean  free  path 
and  not  the  capillary  geometry.  For  each  of  these  five  cases  we 
consider  two  different  wall  temperature  distributions:  linear  and 
step-wise.  In  the  second  and  third  sets  of  calculations,  the  influ¬ 
ence  of  different  wall  temperature  distributions  is  examined  for 
a  channel  and  a  tube,  respectively,  that  connects  reservoirs  with 
a  temperature  ratio  7//T}/  =  2  and  zero  pressure  drop  (cases  Il-a 
to  III-c). 

DSMC  Method 

The  DSMC  method  [6]  has  been  applied  in  this  work  to 
obtain  numerical  solutions  for  short  microchannel  flows.  The 
DSMC  method  is  a  statistical  approach  for  solution  of  the  Boltz¬ 
mann  equation,  the  governing  equation  of  the  rarefied  flows.  The 
DSMC  method,  initially  proposed  by  G.  A.  Bird  in  the  early 
1960s,  has  become  the  most  powerful  and  accurate  method  of 
numerical  modeling  of  complex  rarefied  gas  flows.  [6]  The  fun¬ 
damental  principle  of  the  DSMC  method  is  the  splitting  of  the 
dynamics  of  molecular  motion  during  a  time  step  At  into  two 
sequential  stages:  free-flight  of  molecules  and  intermolecular  bi¬ 
nary  collisions.  The  DSMC  method  is  non-stationary  in  nature, 
but  can  be  used  to  solve  stationary  problems  as  well.  Implemen¬ 
tation  of  the  DSMC  method  usually  implies  the  discretization 
of  the  flow  domain  into  a  grid  of  cells.  The  size  of  computa¬ 
tional  cells  should  be  sufficiently  small  so  that  the  change  in  gas 
dynamic  properties  across  each  cell  is  small.  In  each  cell,  the 
Knudsen  number  based  on  the  cell  size  should  be  larger  than  one. 
When  the  cell  size  in  a  simulation  is  too  big,  macroscopic  gradi¬ 
ents  are  typically  underpredicted  and  the  solution  corresponds  to 
an  artificially  larger  Knudsen  number.  The  time  step  in  the  sim¬ 
ulation  is  usually  selected  such  that  At  =  min(%i,%res),  where 
is  the  mean  time  between  collisions,  xres  is  the  mean  residence 
time  in  a  cell,  so  that  the  molecules  do  not  cross  more  than  one 
cell  during  a  time  step.  At  each  time  step,  the  boundary  condi¬ 
tions  of  the  flow  problem  are  modeled  through  particle  injection 
at  the  domain  boundaries  and  collisions  with  solid  surfaces.  Af¬ 
ter  steady  flow  is  reached,  sampling  of  macroparameters  within 
each  cell  is  performed  for  a  time  period  long  enough  to  avoid 
statistical  scatter. 

The  DSMC -based  software  SMILE  [7]  is  used  for  all  DSMC 
computations.  The  code  uses  the  majorant  frequency  scheme  [8] 
of  the  DSMC  method  for  modeling  of  collisional  process.  The 
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intermolecular  potential  is  assumed  to  be  the  variable  soft  sphere 
model  [9].  The  Larsen-Borgnakke  model  [10]  with  temperature- 
dependent  Z,  and  Zv  and  discrete  rotational  and  vibrational  en¬ 
ergies  is  used  for  the  energy  exchange  between  the  translational 
and  internal  modes.  At  the  gas-wall  interface  the  the  Maxwell 
model  is  used  with  the  fully  diffuse  accommodation  of  tangen¬ 
tial  momentum  and  energy. 

The  background  cell  size  in  2D  DSMC  calculations  was 
l/40th  of  the  channel  height,  the  the  cell  linear  dimension  is 
equal  or  less  than  the  gas  mean  free  path  in  all  considered  cases. 
The  total  number  of  background  cells  in  each  of  the  considered 
cases  was  96,000.  The  average  number  of  simulated  molecules 
per  background  cell  was  10.  The  total  number  of  simulated 
molecules  in  the  domain  was  about  1  million. 


BGK  Model  Kinetic  Equation  Solution 

The  computational  cost  of  DSMC  simulations  of  low-speed 
flows  in  long  channels  becomes  prohibitive  due  to  low  signal- 
to-noise  ratio  [11].  This  is  why  for  modeling  of  such  flows  it  is 
reasonable  to  consider  the  solution  of  the  steady-state  Boltzmann 
equation  with  the  collision  operator  simplified  using  BGK  model 
[12], 

df  df 

U^-+V^-=v(f0-f)  (1) 

dr  ay 

where  f(x,y,U,V,W)  is  the  distribution  function,  x  and  v  are 
Cartesian  coordinates,  U,V,W  are  velocity  components,  v  is  the 
collision  frequency  and  /o  is  the  local  equilibrium  Maxwell- 
Boltzmann  distribution  function  which  depends  on  the  moments 
of  the  distribution  function  /.  This  is  a  non-linear  integrodiffer- 
ential  equation  that  uses  a  simplified  form  of  the  collision  integral 
compared  to  the  Boltzmann  equation,  and  is  applicable  to  model 
gas  flows  with  arbitrary  Knudsen  numbers  and  degree  of  flow 
nonequilibrium. 

The  BGK  model  kinetic  equation  can  be  solved  using  a 
finite-difference  method  and  standard  discrete-ordinate  method 
procedures  [13, 14]  for  interpolation  of  the  distribution  function 
in  velocity  space.  First,  for  the  two-dimensional  problems,  re¬ 
duced  distribution  functions  are  introduced: 

g  =  J  fdW,  h  =  J  fW2dW 

One  can  introduce  the  polar  coordinates  in  the  velocity  space  as 
follows: 

f/  =  5sin(|),  V  =  5cos(|). 


The  discrete  speeds  and  velocity  angles  (])a: 

n  =  So,  nU  =  ££>8P0S5sin(])0,g8o 

8o  8o 


nV  =  ££P8P0S8cos(])0,g8o,  ^nT  =  ^^^o^So  +  ^Sa) 
8o  z 

Finally,  the  resulting  system  of  equations 

5§Sin(')o“§7'  +5Sc°s(|)a^^  =  v(g8o -g&c) 


c  •  a  ^!8o  .  c  i  ^!8o  /,  0  /  \ 

s5  S1  n  <t>0  +  .S’5  cos  <|>o  =v(/!ga-/z8a) 

is  solved  by  finite-difference  method.  The  Gauss-Hermite  half 
range  quadrature  [15]  of  order  16  and  three-eights  rule  are  used 
for  integration  over  the  speed  and  velocity  angle,  respectively. 

Short  Channel  Flow  under  Temperature  and  Pressure 
Differences 

The  flow  in  the  capillary  section  of  a  Knudsen  compressor 
is  driven  by  the  temperature  gradient  in  the  presence  of  a  small 
pressure  gradient.  In  order  to  study  separately  the  temperature 
effect,  a  temperature  drop  is  assumed  across  the  channel,  with 
the  density  in  the  reservoirs  being  equal.  Such  flows  through 
finite  length  channels  in  the  presence  of  both  temperature  and 
pressure  differences  are  examined  in  this  section.  The  results 
presented  in  this  section  correspond  to  Cases  I-a  to  I-e  in  Table  1 . 
In  all  of  the  considered  cases,  the  pressure  ratio  Pi /  Pn  =  1.5  and 
temperature  ratio  is  7//7}/=1.5.  Since  the  pressure  ratio  in  these 
cases  exceeds  the  thermomolecular  pressure  ratio,  there  is  a  non¬ 
zero  mass  flow  from  the  left  reservoir  to  the  right  one.  For  each 
set  of  reservoir  conditions  two  different  temperature  distributions 
are  simulated,  (1)  a  linear  distribution  from  7}  to  7//,  and  (2)  a 
step-wise  distribution  Tw(x)  =  7}  for  x  <  L/ 2;  Tw(x)  =  7}/  for 
x  >  L/2.  The  stepwise  temperature  distribution  may  occur  in 
Knudsen  compressor  if  the  transpiration  membrane  is  composed 
of  the  materials  with  different  thermal  conductivities. 

The  calculated  gas  translation  temperature  and  X- 
component  of  the  velocity  fields  are  shown  in  Figs.  2  and 
3,  respectively,  for  Cases  I-c  (Kn/7=1)  and  I-e  (Kn//=0.05).  The 
temperature  gradients  are  less  steep  for  Case  I-c  than  for  Case 
I-e  due  to  larger  mean  free  path.  The  end  effects  are  apparent 
in  the  velocity  field  and  result  in  an  effective  elongation  of  the 
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channel.  The  most  important  conclusion  that  may  be  drawn 
from  these  figures  is  that  the  temperature  distribution  does  not 
impact  the  flow  outside  of  the  channels.  This  is  true  for  both 
Knudsen  numbers  considered.  Moreover,  the  velocity  fields  for 
the  linear  and  step-wise  temperature  distributions  are  close  near 
the  surface,  but  differ  significantly  at  the  centerline. 

Comparison  of  density  profiles  along  the  centerline  for  the 
above  cases  is  given  in  Fig.  4.  General  behavior  of  density  pro¬ 
files  is  similar  for  different  Knudsen  numbers  and  temperature 
variation,  with  density  being  lower  in  the  high-temperature  part, 
and  higher  in  the  low-temperature  part.  This  is  related  to  a  rela¬ 
tively  lower  thermal  velocities  of  the  latter  and,  as  the  result,  their 
longer  residence  time  and  larger  densities.  Similarly,  the  popu¬ 
lation  of  molecules  in  the  high  temperature  region  is  affected  by 
short  residence  times  in  that  region.  This  effect  is  obviously  more 
pronounced  in  the  step  distribution  cases.  The  molecular  colli¬ 
sions  increase  the  deviation  of  density  from  its  average  value. 

The  X-component  of  velocity  profile  along  the  center- 
line  is  plotted  in  Fig.  5.  The  pressure-driven  Poiseuille  and 
temperature-driven  thermal  creep  flows  represent  two  compet¬ 
ing  flow  mechanisms  in  these  cases.  Due  to  the  thermal  creep 
flow  from  the  cold  container  to  the  hot  one,  the  X-component  of 
velocity  is  lower  in  Cases  I-e  Clin)  and  I-e  (step)  than  in  the  case 
of  no  temperature  gradient.  In  the  case  of  step-wise  wall  temper¬ 
ature  distribution,  the  resulting  velocity  profile  in  the  channel  is 
non-mono  tonic,  and  has  a  strong  gradient  in  the  step  region  due 
to  thermal  creep. 

The  calculated  mass  flow  rates  for  Cases  I-a  to  I-e  are  listed 
in  Table  2.  The  flow  coefficients  Qp  and  Qj  are  used  here  that 
are  usually  employed  for  characterization  of  the  Poiseuille  and 
thermal  creep  flows  between  parallel  plates  for  the  conditions  of 
small  temperature  and  pressure  gradients.  The  flow  coefficients 
are  defined  as  follows: 


Op  = - .  Or  = - - 

(^2RT0)powh2±^-  {\V2RTQ)p0wh2jr^ 

(2) 

where  M  is  the  mass  flow  rate  (kg/s),  po,  po ,  and  To  are  the  av¬ 
erage  density,  pressure,  and  temperature,  respectively,  and  w,  h, 
L  are  channel  width,  height  and  length,  respectively.  There  is 
no  significant  difference  in  the  mass  flow  rate  between  the  linear 
and  stepwise  temperature  distributions.  In  fact,  the  difference  in 
mass  flow  rate  is  less  than  2  per  cent  for  all  Knudsen  numbers. 


Influence  of  Wall  Temperature  Distribution 

Temperature  gradient  driven  flows  are  very  sensitive  to  ther¬ 
mal  boundary  conditions.  For  practical  microdevices  it  is  gener¬ 
ally  difficult  to  attain  idealized  thermal  boundary  requirements. 
Thus,  numerical  methods  for  reliably  estimating  the  magnitude 


of  the  consequences  induced  by  non-ideal  conditions  are  impor¬ 
tant. 

As  an  illustration  of  the  usefulness  of  this  approach  consider 
the  recent  development  by  Young  [16]  of  a  prototype  Knudsen 
Compressor.  Radiant  energy  was  used  to  generate  temperature 
gradients  in  the  carbon  doped  aerogel  membranes  that  are  the 
active  elements  in  the  compressor’s  pumping  stages.  Due  to  fi¬ 
nite  photon  mean  free  paths  in  the  carbon  doped  aerogel  mem¬ 
branes,  the  wall  temperatures  in  the  membrane’s  flow  passages 
had  a  maximum  near  to,  but  not  at,  the  exit  surface  of  the  mem¬ 
brane.  This  is  the  motivation  for  the  calculations  in  this  section 
with  the  non-monotonic  wall  temperature  variation  of  case  III-c 
(see  table  1)  constructed  to  mimic  the  essentials  of  a  membrane 
in  Young’s  experiments. 

The  influence  of  wall  temperature  distribution  on  the  ther¬ 
mal  transpiration  flow  in  short  microchannel  was  considered  for 
a  channel  with  length-to-height  ratio  of  5  at  Knu= 0.2.  To  study 
this  effect,  the  pressure  was  assumed  equal  in  the  two  reservoirs, 
and  the  temperature  drop  was  set  to  300  K.  The  three  cases  of 
temperature  distribution  are  (i)  the  linear  temperature  variation 
along  the  channel  wall  (case  ITa),  (ii)  a  step-wise  temperature 
variation  (case  Il-b)  and  (iii)  a  non-monotonic  distribution  typ¬ 
ical  for  a  radiatively  heated  nanomembrane  Knudsen  compres¬ 
sor  [16]  (case  II-c).  A  non-monotonic  wall  temperature  profile 
has  a  maximum  of  750  K  at  one-fourth  of  the  channel  length. 
The  three  wall  temperature  distributions  are  plotted  in  Fig.  6. 

The  change  in  the  wall  temperature  distribution  significantly 
impacts  the  flow  velocity  (see  Fig.  8).  Note  that  only  a  portion 
of  computational  domain  is  shown  to  illustrate  the  details  of  the 
flow  inside  the  channel.  This  impact  is  noticeably  stronger  that 
the  one  shown  in  the  previous  section,  and  is  observed  both  at  the 
centerline  and  near  the  surface. 

The  influence  of  wall  temperature  distributions  inside  the 
channel  on  the  temperature  field  is  illustrated  in  Fig.  7.  The  gas 
temperature  generally  follows  the  surface  temperature  inside  the 
channel,  with  the  end  effects  being  visible  even  at  a  distance  from 
the  inlet/outlet  equal  to  the  channel  length.  The  end  effects  are 
somewhat  smaller  for  the  stepwise  distribution,  and  largest  for 
the  non-monotonic  profile. 

Figure  9  shows  the  calculated  pressure  profiles  along  the 
centerline  for  the  above  temperature  distributions.  The  qualita¬ 
tive  behavior  is  similar  for  the  linear  and  stepwise  distribution, 
with  pressure  larger  in  the  hot  part,  and  smaller  in  the  cold  part. 
The  non-monotonic  case  significantly  differs  from  the  first  two, 
since  there  is  no  significant  pressure  minimum  in  the  cold  part. 
Also,  there  is  a  pressure  increase  inside  the  channel  of  about  10 
percent  for  the  monotonic  profile  due  to  the  reverse  transpiration. 

The  calculated  mass  flow  rate  and  flow  coefficient  Qp  for 
the  three  cases  are  listed  in  Table  3.  The  mass  flow  rate  in  case  of 
stepwise  temperature  variation  is  only  slightly  lower  (by  ss  3%) 
than  for  the  linear  temperature  variation.  However,  the  reverse 
transpiration  in  the  non-monotonic  temperature  distribution  case 
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results  in  significant  degradation  of  mass  throughput  (by  «  1 8%). 

Similar  influence  of  the  wall  temperature  distribution  is  ob¬ 
served  for  a  cyllindrical  tube  case  (see  Table  4).  The  decrease  of 
the  mass  flow  rate  is  close  to  that  for  a  two-dimensional  case  and 
is  equal  to  about  2%  and  17%  for  stepwise  and  non-monotonic 
wall  temperature  distributions,  respectively. 

Validation  of  the  BGK  Solver 

Results  of  the  validation  study  of  a  channel  flow  modeling 
using  a  BGK  model  kinetic  equation  are  presented  in  this  section. 
The  solutions  obtained  using  the  BGK  solver  are  compared  with 
the  solution  of  Navier-Stokes  equations  for  long  channels  at  low 
Knudsen  number  and  a  DSMC  solution  for  a  short  channel  in  the 
transitional  regime. 

1)  Comparison  with  Navier-Stokes  solution  The  flow  of  ni¬ 
trogen  in  a  channel  with  a  length-to-height  ratio  of  30  was  cal¬ 
culated  for  a  pressure  ratio  of  1 .5  and  a  constant  temperature  of 
the  wall  of  300  K.  The  temperature  in  the  upstream  and  down¬ 
stream  reservoirs  is  equal  to  300  K.  The  Knudsen  number  based 
on  the  channel  height  and  the  mean  free  path  in  the  downstream 
reservoir  is  0.05. 

The  X-component  of  velocity  field  is  shown  in  Fig.  10  ob¬ 
tained  using  three  different  uniform  spatial  grids  in  the  BGK 
solver.  The  grid  dimensions  are:  (i)  101x51,  (ii)  51x26  and 
(iii)  26x14.  The  analytical  approximate  solution  [17]  of  Navier- 
Stokes  equations  for  a  long  channel  flow  with  slip  boundary  con¬ 
ditions  is  also  plotted  in  Fig.  10.  The  discrete  ordinate  solution  of 
BGK  model  kinetic  equation  with  coarser  grids  results  in  a  lower 
velocity  in  the  whole  computational  domain.  Note  that  a  simi¬ 
lar  dependence  of  results  on  the  number  of  simulated  particles 
and  cells  in  the  DSMC  method  was  observed  [18],  The  solution 
obtained  for  a  101x51  grid  is  in  excellent  agreement  with  the 
Navier-Stokes  solution. 

2)  Comparison  with  a  DSMC  solution  The  flow  of  nitrogen 
in  a  short  channel  with  a  length-to-height  ratio  of  5  was  calcu¬ 
lated  for  a  pressure  ratio  of  1.5  and  a  constant  wall  temperature 
of  300  K.  The  Knudsen  number  based  on  the  channel  height  and 
the  mean  free  path  in  the  downstream  reservoir  is  0.2,  resulting 
in  a  transitional  flow  regime  in  the  channel. 

Figure  1 1  shows  a  comparison  of  X-component  of  velocity 
profiles  in  the  channel  obtained  by  the  BGK  solver  (top)  and  a 
DSMC  solver  SMILE  (bottom).  Overall,  the  two  solutions  are 
in  a  very  good  agreement  with  the  maximum  difference  in  the 
velocity  values  at  the  axis  of  about  2%  which  is  within  the  error 
of  both  numerical  methods. 

Thus,  the  area  of  applicability  of  the  BGK  model  kinetic 
equation  for  numerical  modeling  of  microchannel  flows  spans  a 
wide  range  of  flow  regimes:  from  near-continuum  to  rarefied, 
as  well  as  different  geometries:  both  small  and  large  length-to- 
height  ratios,  and  arbitrary  pressure  gradients  between  the  reser¬ 
voirs. 


Conclusions 

The  direct  simulation  Monte  Carlo  method  and  the  discrete 
ordinate  method  for  the  solution  of  the  BGK  model  were  used 
to  numerically  study  rarefied  gas  flows  through  two-dimensional 
finite  length  channels.  A  wide  range  of  flow  rarefaction  regimes 
was  examined,  from  Kn=0.05  to  50,  in  the  presence  of  pressure 
and/or  temperature  gradients.  Three  different  wall  temperature 
distributions  were  considered,  namely,  linear,  step-wise,  and  a 
non-monotonic  profile  typical  for  a  radiantly  heated  Knudsen 
compressor  membrane. 

The  influence  of  wall  temperature  distribution  on  the  ther¬ 
mal  transpiration  flow  in  a  short  channel  at  Knu= 0-2  was  consid¬ 
ered.  The  non-monotonic  wall  temperature  distribution  case  sig¬ 
nificantly  differs  from  the  first  two.  There  is  a  pressure  increase 
inside  the  channel  of  about  10  percent  for  the  non-monotonic 
profile  due  to  the  reverse  transpiration.  The  mass  flow  rate  in  case 
of  stepwise  temperature  variation  is  only  slightly  lower  (by  3%) 
than  for  the  linear  temperature  variation.  However,  the  reverse 
transpiration  caused  by  the  non-monotonic  wall  temperature  dis¬ 
tribution  results  in  significant  degradation  of  mass  flow,  by  about 
18%.  Similar  influence  of  the  wall  temperature  distribution  is 
observed  for  a  cyllindrical  tube  case. 

The  application  of  the  numerical  techniques  discussed  in  this 
paper  to  detailed  investigations  and  optimizations  of  microscale 
thermal  creep  flows  will  be  extremely  valuable.  The  work  re¬ 
ported  here  is  a  first  look  at  describing  practical  thermal  gradient 
driven  or  thermal  creep  flows.  While  the  available  results  are  at 
present  both  scarce  and  preliminary,  they  do  validate  in  this  case 
what  was  anticipated,  and  provide  useful  quantitative  estimates 
of  the  consequences  of  the  particular  wall  temperature  maximum 
chosen  for  the  calculations. 
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Table  1 .  Summary  of  cases  considered 
Flow  Conditions  Designation 


Pressure  and  Temperature-Driven  Flow 


L/h=5,  N2,  Tr- 

=600  K,  T//=300  K,  ^-=1.5 

pii 

Kn, /=50 

case  I-a 

Knu= 3 

case  I-b 

Knu=  1 

case  I-c 

Kn,  i=Q.2 

case  I-d 

Knu= 0.05 

case  I-e 

Effect  of  Wall  Temperature  Distribution 
Two-dimensional  channel: 

L/h=5,  Kn//=0.2,  N2,  7)=600  K,  7)/=300  K,  -^=1.0 
linear  case  Il(lin) 

step  case  II(step ) 

non-monotonic  case  Il(nmon) 

Axisymmetric  tube: 

L/a=5,  Kn//=0.2,  N2,  7}=600  K,  7}/=300  K,  ^=1.0 
linear  case  Ill(lin) 

step  case  III( step ) 

non-monotonic  case  III( nmon ) 

Table  2.  Calculated  mass  flow  rate  and  flow  coefficients  for  pressure  and 
temperature  driven  flow  in  short  channel. 


Case 

m 

(kg/s) 

Qp 

Tp  x  107 
( kg/s-Pa) 

I-a  (lin) 

0.119xl0~7 

0.533 

2.17 

I-a  (step) 

0.119xl0~7 

0.533 

2.17 

I-b  (lin) 

0.229xl0~6 

0.615 

2.5 

I-b  (step) 

0.229x  10~6 

0.615 

2.5 

I-c  (lin) 

0.812x  10~6 

0.727 

2.96 

I-c  (step) 

0.824x  10~6 

0.738 

3.0 

I-d  (lin ) 

0.672x  10~5 

1.20 

4.90 

I-d  (step) 

0.677xl0~5 

1.21 

4.93 

I-e  (lin) 

0.555  xl0~4 

2.49 

10.11 

I-e  (step) 

0.554x  10~4 

2.48 

10.09 
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Table  3.  Calculated  mass  flow  rate  and  flow  coefficient  for  temperature- 
driven  flow  in  short  channel,  Case  II. 


linear 


Temperature  Distribution 

m  (kg/s) 

Qt 

linear 

1.580E-6 

0.1189 

stepwise 

1.532E-6 

0.1153 

non-monotonic 

1.292E-6 

0.0972 

Table  4.  Calculated  mass  flow  rate  and  flow  coefficient  for  temperature- 
driven  flow  in  short  tube,  Case  III. 


Temperature  Distribution  lit  (kg/s)  Qj 

linear  1.06E-9  0.1269 

stepwise  1.04E-9  0.1245 


Kn-l 


step 


mumming 


■wwj  j 

linear 

Kn-0.05 

step 

T,  K 

300.00 

307.50 
315.00 

322.50 
330.00 

337.50 
345.00 

352.50 
360.00 

367.50 
375.00 

382.50 
390.00 

397.50 
405.00 

412.50 
420.00 

427.50 

1435.00 
442.50 
450.00 


Figure  2.  Translational  temperature  field  for  Cases  l-c  and  l-e. 


non-monotonic  8.8E-10  0.1054 


Figure  3.  X-component  of  velocity  field  for  Cases  l-d  and  l-e. 
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Figure  4.  Normalized  density  profile  along  the  channel  for  different  tem¬ 
perature  variations  and  Knudsen  numbers  (Cases  l-c  and  l-e). 


Figure  6.  Wall  temperature  distributions  for  Case  II. 


Figure  5.  X-component  of  velocity  profile  along  the  channel  for  different 
temperature  distributions  and  Knudsen  numbers  (Cases  l-c  and  l-e). 


Figure  7.  Translational  temperature  field,  Kn//=0.2,  L/h=5  for  different 
temperature  variations,  Case  II 
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lineal • 


step 


non¬ 

monotonic 


U,  m/s 

-10. 

-9.5 

-9.0 

-8.5 

-8.0 

-7.5 

-7.0 

-6.5 

-6.0 

-5.5 

-5.0 

-4.5 

-4.0 

-3.5 

-3.0 

-2.5 

-2.0 

-1.5 

-1.0 


Figure  8.  X-component  of  velocity  field,  Kn//=0.2,  L/h=5  for  different 
temperature  variations,  Case  II 
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U,  m/s 

I  18.0 

1  17.0 
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-  9.0 

-  6.0 

-  4.0 

J3.0 
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Figure  10.  Comparision  of  X-component  of  velocity  field  calculated  us¬ 
ing  Navier-Stokes  equations  (top)  and  BGK  model  with  different  meshes, 
Kn//=0.05,  L/h=30. 


Figure  9.  Pressure  profile  along  the  channel,  Kn//=0.2,  L/h=5  for  differ¬ 
ent  temperature  variations,  Case  II.  Figure  1 1 .  Comparison  of  X-component  of  velocity  field  calculated  using 

BGK(top)  and  DSMC  (bottom)  models,  Kn//=0.2,  L/h=5. 


9 


Approved  for  public  release;  distribution  unlimited.  Copyright  ©  2005  by  ASME 


